Morphology transitions induced by chemotherapy in carcinomas "in situ 
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Recently, we have proposed a nutrient-limited model for the avascular growth of tumors including 
cell proliferation, motility and death [Q, that, qualitatively reproduces commonly observed mor- 
phologies for carcinomas in situ. In the present work, we analyze the effects of distinct chemothera- 
peutic strategies on the patterns, scaling and growth laws obtained for such nutrient-limited model. 
Two kinds of chemotherapeutic strategies were considered, namely, those that kill cancer cells 
and those that block cell mitosis but allows the cell to survive for some time. Depending on the 
chemotherapeutic schedule used, the tumors are completely eliminated, reach a stationary size or 
grow following power laws. The model suggests that the scaling properties of the tumors are not 
affected by the mild cytotoxic treatments, although a reduction in growth rates and an increase in 
invasiveness are observed. For the strategies based on antimitotic drugs a morphological transition 
in which compact tumors become more fractal under aggressive treatments was seen. 



PACS numbers: 87.19+e, 87.18. Hf, 87.18.Bb, 87.15.Vv 



I. INTRODUCTION 

Cancer is the uncontrolled cellular growth in which 
neoplastic cells invade adjacent normal tissues and give 
rise to secondary tumors (metastasis) on tissues or or- 
gans distant from its primary origin |y. A cancer that 
remains confined within a messed surface, without break 
of the underlying basement membrane, is referred to as 
carcinoma "in situ" . "In situ" carcinoma is characterized 
by intense cytological atypia, necrosis, and frequent and 
abnormal mitosis, the tumor cells being arranged in vari- 
ous distinctive patterns || . A malignant tumor is derived 
from a single cell that, years before the tumor becomes 
clinically detectable, began an inappropriate pathway of 
proliferation Although cancers are extremely diverse 
and heterogeneous, a small number of pivotal steps as- 
sociated to both deregulated cell proliferation and sup- 
pressed cell death are required for the development of 
any and all tumors. Indeed, all neoplasms evolve ac- 
cordingly to an universal scheme jj^, g|. In the struggle 
against cancer, surgical removal, chemotherapy and/or 
radiotherapy are the most commonly used treatment for 
the complete eradication of the tumor mass. Nowadays, 
new approaches, such as immunotoxins Q|, gene ||, anti- 
angiogenic H and virus |]Io|| therapies, are being devel- 
oped and have been successfully used in the treatment of 
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several kinds of experimental and human tumors. 

Mathematical models are always used as a tentative 
for describing cancer growth. In particular, numerous 
models based on classical reaction-diffusion equations 
have been proposed to investigate the growth of tumor 
spheroids [ pd] , cancer progress and its interaction with 
the immune system |12], and the tumor angiogenesis 
jD| O. Fractal growth patter ns in gliomas (brain tu- 
mors) were recently investigated by Sander and Deis- 
boeck Jl5| . Scalerandi et. ai, using the local interaction 
simulation approach (LISA), formulated models for the 
evolution of avascular tumors |l6[ ] and tumor cords pj] 
under nutrient competition. They also analyzed the ef- 
fects of tumor vascularization [H on cancer growth. Re- 
cently, we have studied a diffusion-limited model for the 
growth of carcinomas "in situ" in which cell proliferation, 
motility and death are locally controlled by growth fac- 
tors ^(|. This model generates compact, connected 
and disconnected morphologies characterized by Gom- 
pertzian growth in time and distinct scaling laws for the 
tumor interfaces. In order to generate papillary and ram- 
ified morphologies found in many epithelial cancers and 
trichoblastomas, we were led to analyze the effects of nu- 
trient competition in cancer development |Q. 

In addition to the vast literature dedicated to tu- 
mor growth modeling, many research papers addressing 
cancer therapies have recently been published. Indeed, 
cancer cell kinetics under treatments using antimitotic 
agents |2^| , radiotherapy (23), virus that replicate 
selectively in tumor cells J24J, anti-angiogenic chemicals 
PH , as wells as the effects of tumor drug resistance and 
tumor vasculature on chemotherapies p6[, were studied 
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using mathematical models. 

In this paper it is analyzed the effects of distinct 
chemotherapeutic strategies in the model for the growth 
of avascular tumors proposed by us JlJ . Specifically, we 
are interested in possible changes in the tumor patterns, 
scaling and growth laws reported on the original model 
triggered by chemotherapies. The paper is organized as 
follows. In section [n] the nutrient-limited model for can- 
cer growth and its main results are reviewed. In section 



III, chemotherapeutic treatments that aim to kill can- 
cer cells are introduced in the framework of the present 
model and their effects discussed. In section [Fv|, a model 
for chemotherapy with antimitotic agents is considered. 
Finally, we draw some conclusions in section [v| 



II. THE NUTRIENT-LIMITED MODEL 

Our nutrient-limited model combines macroscopic 
reaction-diffusion equations, describing the nutrient field 
concentration, with microscopic stochastic rules govern- 
ing the actions of individual tumor cells. The basic prin- 
ciples included in the model are cell proliferation, motil- 
ity and death as well as competition for nutrients among 
cancer and normal cells. Also, the nutrient concentration 
field locally controls cell division, migration and death. 



A. The model 

1. The tissue 

The tissue, represented by a square lattice of size 
(L + 1) x (L + 1) and lattice constant A, is fed by a 
single capillary vessel at x = 0, the top of the lattice. 
The capillary is the unique source from which nutrients 
diffuse through the tissue towards the individual cells. 
Although a tumor mass is composed of different cell sub- 
populations ||, the model considers only three types: 
normal, cancer and tumor dead cells. Any site, with co- 
ordinates x = (iA,jA), i,j = 0, 1,2, ...,L, is occupied 
by only one of these cell types. In contrast to the normal 
cells, one or more cancer cells can pile up in a given site. 
In turn, dead tumor cells are inert. Thus, each lattice site 
can be thought of as a group of cells in which the normal, 
dead and cancer cell populations assume one of the pos- 
sible values a n (x) = (Jd{x) = 0, 1 and a c (x) = 0, 1, 2, . . ., 
respectively. Accordingly the theory of the monoclonal 
origin of cancer [Q, a single cancer cell at y = LA/2 and 
at a distance X from the capillary is introduced in the 
normal tissue. Periodic boundary conditions along the 
horizontal axis are used. The row i — represents the 
capillary vessel and the sites with i = L + 1 constitute 
the external border of the tissue. This geometry is par- 
ticularly adequate to describe the growth of carcinomas 
(epithelial tumors) "in situ" because the present model 
consider that the tumor mass receive nutrients only by 
diffusion from the capillary vessel. However, the model 



can be extended to others cancers. 



2. The nutrients 

The nutrients are divided into two classes: essential to 
cell proliferation, those necessary to DNA synthesis, and 
nonessential to cell division. The essential and nonessen- 
tial nutrients are described by the concentration fields 
N(x,t) and M(x,t), respectively. These nutrient fields 
obey the dimensionless diffusion equations (see ref. Q| 
for the complete variable transformations) 



dN 
~dt 



y^N - a 2 Na n - X N a 2 Na, 



and 



dM 9 , 
at 



Am a Ma c 



(1) 



(2) 



in which differentiated nutrient consumption rates for 
normal and cancer cells by factors Xn and Am are as- 
sumed. It is important to notice that the model admits 
the simplest form for nutrient diffusion, i. e., linear equa- 
tions with constant coefficients. Also, Ajv > Am is used, 
reflecting the larger affinity of cancer cells for essential 
nutrients. 

The boundary conditions satisfied by the normalized 
nutrient concentration fields are N(x — 0) = M(x — 
0) = 1, representing the continuous and fixed supply of 
nutrients provided by the capillary vessel; N(y = 0) = 
N(y = L) and M(y = 0) = M(y — L), corresponding 
to the periodic boundary conditions along the y-axis; at 
last, Neumann boundary conditions, dN(x — L)/dx = 
dM(x = L)/dx = 0, are imposed to the border of the 
tissue. 



3. Cell dynamics 

In our original model, each tumor cell, randomly se- 
lected with equal probability, can carry out one of three 
actions: division, migration or death. However, in the 
present work we consider just the accommodation that 
happens during cell mitosis among the daughter cells. 
Consequently, each tumor cell can carry out one of two 
actions. 

1. Division. Cancer cells divide by mitosis with prob- 
ability Pdiv If the chosen cell is inside the tu- 
mor, its daughter will pile up at that site, and 
<j c (x) — * a c (x) + 1. Otherwise, if the selected cell is 
on the tumor border, its daughter cell will occupy 
at random one of their nearest neighbor sites x' con- 
taining a normal or a necrotic cell and, therefore, 
a c (x') = 1 and o- n ,d{x') — 0. The mitotic proba- 
bility P c u v is determined by the concentration per 
cancer cell of the essential nutrients N present on 
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the microenvironment of the selected cell: 

2" 



Pdiv{x) = 1 - exp 



N(x) 



a c {x) 9 dv 



(3) 



The Gaussian term is included in order to produce a 
sigmoid curve saturated to the unity, and the model 
parameter Odiv controls the shape of this sigmoid. 

2. Death. Cancer cells die with probability Pdei- 
Thus, a c (x) — > a c (x) — 1 and ad(x) — 1 when a c 
vanishes. The cell death probability P del is deter- 
mined by the concentration per cancer cells of the 
nonessential nutrients M present on the microenvi- 
ronment of the selected cell: 



Pdei(x) = exp 



M(x) 

<Jc(x) 



<l(d 



(4) 



i. e., a Gaussian distribution whose variance de- 
pends on the model parameter 6 del- 

The biological basis for these cell dynamic rules can 
be found in reference Q. However, it is worthwhile to 
notice that from the point of view of the so-called kinetic 
cellular theory, which provides a general framework for 
the statistical description of the population dynamics of 
interacting cells 12 , the local probabilities Pdiv and Pdei 
can be thought as an effective kinetic cellular model. 



normal and cancer cells, controlled by the model parame- 
ters a, Xn and Am, plays a central role in morphology se- 
lection. For small values of these parameters, correspond- 
ing to growth conditions in which individual cells demand 
small nutrient supplies, the patterns tend to be compact 
and circular (Fig. ^a). However, if the mitotic rate of 
cancer cells is small due to the large amount of nutrients 
demanded for cell division, generating a significant com- 
petition for nutrients, these compact patterns progres- 
sively assume papillary- like morphologies (Fig. gb). At 
high nutrient consumption rates these papillary patterns 
become the rule and, for low cancer cell division, con- 
tinuously transform into thin tips, filaments or chords of 
cells (Fig. ||c). Also, the model generates patterns with 
a inner core of died cells for high nutrient consumption 
or cell division rates (Fig. ||d). As observed in "in vivo" 
tumors and "in vitro" multicell spheroids Jll|] , these sim- 
ulated patterns consist of three distinct regions: a central 
necrotic core, an inner rim of quiescent cancer cells and 
a narrow outer shell of proliferating cells. 

The tumor patterns generated by our nutrient-limited 
model were characterized by its gyration radius R g , total 
number of cancer cells Nc, and number of sites on tumor 
periphery S (including the surface of holes, if any). The 
gyration radius R g is defined as 



1/2 



(5) 



4- Computer implementation 

The growth model simulations were implemented us- 
ing the following procedure. At each time step T, Eqs. 
(Q) and (|^) are numerically solved in the stationary 
state (dN/dt = dM/dt — 0) through relaxation meth- 
ods. Provided the nutrient concentration at any lattice 
site, Nc(T) cancer cells are sequentially selected at ran- 
dom with equal probability. For each one of them, a 
tentative action (division or death) is randomly chosen 
with equal probability and the time is incremented by 
AT = 1/N C (T). The selected cell action will be im- 
plemented or not according to the correspondent local 
probabilities determined by Eqs. (||) or (jl]). At the end 
of this sequence of Nc(T) tentatives, a new time step 
begins and the entire procedure (solution of the diffusion 
equations and application of the cell dynamics) is iter- 
ated. The simulations stop if any tumor cell reaches the 
capillary vessel. 



B. Main results 

The model reproduces commonly observed tumor mor- 
phologies including the papillary, compact and ramified 
patterns shown in Fig. 0. The nutrient consumption by 



where n is the number of sites occupied by the pat- 
tern (cancer or necrotic cells) and r» is the distance of 
the occupied site i from the tumor mass center. These 
quantities could be related to clinically important cri- 
teria such as progress curves, growth rates (volumetric 
doubling time) at given radii, proliferative and necrotic 
fractions of the tumor. 

The progress in time of cancer cell populations for all 
the simulated patterns follows Gompertz curves, 



N C (T) = A exp [- exp(-k(T - T c ))] , (6) 



representing an initial exponential growth that is satu- 
rated in larger times. 

Alternatively, as a function of the number of sites oc- 
cupied by the pattern n, both R g and S obey power law 
scalings given by R g ~ n v and S ~ n a , respectively. For 
solid patterns these exponents are v ~ 0.5 and a ~ 0.5, 
corresponding to effective circular and non fractal pat- 
terns. As the nutrient consumption increases, the pat- 
terns tend to papillary-like shapes for which the exponent 
a increases towards the value 1 and the exponent v varies 
in the range [0.50,0.60], indicating a fractal morphology 
for the tumor. 
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III. CYTOTOXIC CHEMOTHERAPEUTIC 
DRUGS 

The primary aim of the antitumoral treatment with 
chemotherapeutic drugs is to kill or at least to stop the 
proliferation of the cancer cells. In general, the drugs 
should act only on proliferating cells, mainly the cancer 
ones. However, drugs also destroy proliferating normal 
cells promoting several collateral effects ^7j. Indeed, ep- 
ithelial cells from the respiratory and gastro-intestinal 
systems, which frequently reproduce in order to substi- 
tute their dead counterparts, are strongly affected. In 
this section, we analyze a simple chemotherapeutic model 
in which the complex details of the cell-cycling responses 
to the drugs are taking into account by an effective ki- 
netic model. 



The model 





X 










\ Exponential 








\ decay 




Treatment 








beginnig 









FIG. 1: Temporal profile of drugs concentration at the capil- 
lary vessel. 



where Z = 0,1,2... is the number of applied doses, T x is 
the characteristic time for the drug elimination by the or- 
ganism and To is the time at which the treatment begins. 
The functional form of Q(T) is shown in Fig. [I. 

As the nutrients, drugs diffuse from the capillary ves- 
sel toward the individual cells and their concentration 
Q(x,t) at any lattice site at each time step is given by 
the stationary solution of the diffusion equation 



dQ 
dt 



-£ = DqV 2 Q - T 2 Q - \qT 2 Q(j c 



(8) 



The first, second and third terms on the right hand side 
represent drug diffusion, natural degradation and absorp- 
tion by cancer cells, respectively. Again, this equation is 
the simplest one describing the diffusion phenomenon. 
Eq. (||) is written using the same dimensionless variables 
of |l]) and (||). Also, the parameters Qo and Dq can 
be make equal to the unity without generality loss. The 
boundary conditions are the same used for the nutrient 
fields, except at the capillary vessel where the concentra- 
tion at each time step is given by Eq. (0). 

Finally, a single change is introduced into the cell dy- 
namics described in section [S[ an additional chance of 
death occurs whenever a cancer cell enters mitosis. Thus, 
every time a cancer cell divides, each one of the generated 
cells can die with probability 



i , i2 ) (2) = l-exp 



Q(x) 



(TAX) 



n{Q) 
'del . 



(9) 



The parameter 6^ controls the cell sensitivity to the 
drugs. Also, since in the original nutrient-limited model 
the normal cells do not take part on the cell dynamics, i. 
e., they do not divide or die, we disregard the chemother- 
apeutic effects on normal cells. 



As used in previous works 128|, the chemotherapy 
is modelled by a periodic delivery of cytotoxic drugs 
through the same capillary vessel supplying the nutrients 
to the tissue. Several cytotoxic drugs and their properties 
were exhaustively studied such as amsacrine, cisplatin, 
cyclophosphamide, cytarabine, mustine, and anthracy- 
cline [p9| . Here, the numerous barriers involved in tumor 
drug delivery p0| were not considered, and the treatment 
begins when the tumor mass contains Nq cancer cells. 
When a dose is applied, the drug concentration in the 
capillary assumes a value Qo- By hypothesis, this con- 
centration level decays exponentially in time, simulating 
the gradual drug elimination by the organism. New doses 
are periodically administered at time intervals r. So, the 
drug concentration at the capillary, Q(T), is a function 
of time written as 

f 0, if T < T 
Q(T) = { Qoexp[-{T-l T )/T x }, if T + lr<T< , 
[ T + (I + 1)t 

(7) 



B. Results 

The main aim of this work was to investigate the ef- 
fects of treatments on the various morphologies, scaling 
and growth laws observed in the original model. Thus, 
for each morphology, determined by the fixed parameters 
reported in table pL the treatment parameters, namely, r 

(O) 

(the dose period), 9^{ (drug efficiency) and N (tumor 
size at the treatment beginning) were varied. Such pa- 
rameters can be directly tested in the laboratory. The 
remaining parameters (Xq, T and T x ) associated to drug 
diffusion and elimination were also fixed in all the simu- 
lations. 

In Fig. ||, the corresponding patterns for compact, 
papillary, and ramified morphologies with and without 
treatment are shown. In this figure the treatment is not 
able to cease tumor development. As can be seen, the 
morphological tumor patterns do not change under mild 
treatment. However, the regions occupied by the tumors 
are larger than those without treatment. These results 
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Morphology 


Ajv 


Am 


Xq a T 8 div 8 de i 


T x 


Compact 


25 


10 


10 2/L 2/L 0.3 0.03 


4 


Papillary 


200 


10 


10 2/L 2/L 0.3 0.03 


4 


Ramified 


200 


10 


10 3/L 2/L 0.3 0.01 


4 


Necrotic 


50 


25 


10 3/L 2/L 0.3 0.03 


4 



TABLE I: Fixed parameters used in tumor growth simulations 
under chemotherapeutic treatment for each morphology type. 



suggest that the direct attack to the tumor might be an 
inadequate treatment strategy. Indeed, the more invasive 
is the tumor, higher is the possibility of the cancer cells to 
reach the capillary vessel and, consequently, metastasize 
successfully Such result is in agreement with the Israel's 
claim J3l| that cancer cells trigger adaptation mecha- 
nisms in stress circumstances similar to those observed 
in bacterial colonies. Moreover, in our model neither ge- 
netic or epigenetic changes are necessary to explain the 
increase of tumor aggressiveness. It naturally emerges 
from the growth rules. A tumor submitted to succes- 
sive chemotherapeutic treatments which do not lead to 
its complete eradication may progressively become more 
resistant, aggressive and malignant. 

Very similar scaling laws R g ~ n" and S ~ n a for 
the treated tumors were observed {y and a values are 
reported in The small differences in the exponents 
values for the number of tumor peripheral sites vanishes 
at the asymptotical limit of tumor size. This exponents 
invariance suggests that the fractal morphology is a ro- 
bust aspect of these tumors and cannot be changed by 
small perturbations in the cell microenvironment. Obvi- 
ously, the scaling laws for R g and S are not defined for 
the tumors that cease their growth. 

We have also studied the influence of the treatment 
parameters (t, fljS , No) on the tumor growth curves. In 
Fig. H the evolution in time of the total number of cancer 
cells for compact tumors is shown. Depending on dose 
period r the tumors may disappear, saturate their growth 
or progress according to power laws. Actually, the time 
interval between two consecutive doses is a fundamental 
clinical feature determining the treatment success. The 
power law regime observed when r is large means a slow 
progress, neatly contrasting with the exponential growth 
present in the Gompertz law describing tumor progress 
without chemotherapy. The tumor size at the beginning 
of the treatment is another important factor. Indeed, 
the exponents of the power laws characterizing the tumor 
growth are smaller when the initial tumor size is bigger. 
This scenario seems to be universal. Sometimes, as shown 
Fig. |^, compact tumors which began to receive drug 
doses at regular intervals of r = 5 when it contained No = 
10 4 cancer cells were eliminated slower than another one 
that began to be treated with No = 5 X 10 4 cancer cells. 
However, this is not the rule. For example, simulations of 
papillary tumors indicate that smaller tumors are faster 
eliminated when r = 5. Finally, since the growth law 
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(d) 



FIG. 2: Tumor growth patterns generated by the limited- 
nutrient model under chemotherapeutic treatment. The pat- 
terns are drawn in a gray scale where the darker gray regions 
represent higher cancer cell populations, the black points rep- 
resent the sites occupied by necrotic cells, and the white re- 
gions represent the normal tissue. The tissue size is 500 x 500, 
and the initial "cancer seed" is 300 sites distant from the cap- 
illary. The total number of cancer cells depends on tumor 
morphology and reach up to 2 x 10 5 for compact patterns. 
Two typical patterns, without (left) and with soft treatment 
(right) are shown, for (a) compact, (b) papillary, (c) ramified 
morphologies and (d) patterns with a necrotic core. The fixed 
parameters used for generate these morphologies are listed in 
table Mild treatment means that the period between two 
doses is large (r = 20) and, therefore, the tumor grows contin- 
uously. The other parameters are 9^ =0.1 and iVo = 10 4 . 



exponents depend on the parameter sets used, they are 
not universal. Concerning the 9^$ parameter, it just 
modifies the treatment efficiency. 

Fractal tumors are more resistant to treatments. In 
Fig. the growth curves Nq x T are drawn for compact, 
papillary, ramified and necrotic morphologies shown in 
Fig. ^. One can see that the more fractal is the tumor 
larger is the required time to eliminate it. Indeed, the 
lower growth rates of fractal tumors imply in a large frac- 
tion of cancer cells maintained at a quiescent state. So, 
since chemotherapeutic drugs considered in this model 
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IV. ANTIMITOTIC CHEMOTHERAPEUTIC 
DRUGS 



FIG. 3: Growth curves for compact patterns (dashed lines). 
Three dose intervals were tested (r = 5, 10, 20) in tumors with 
two distinct initial sizes: No = 10 4 and A = 5 x 10 4 (inset). 
The exponents of the power laws (slopes of the straight lines) 
are 1.3 (0.98) for the tumors with N = 10 4 (N = 5 x 10 4 ) 
when r = 20, and 0.05 (-0.09) when r = 10. 



only act in the dividing cells, the treatment becomes in- 
efficient when the major fraction of the cancer cells are 
quiescent. 




FIG. 4: Growth curves for distinct tumor morphologies. The 
doses were applied at intervals r = 5 in tumors with initial 
sizes iVo = 10 4 . 

In addition, we have studied tumor patterns exhibit- 
ing a central necrotic core. The results for cancer growth 
and the correspondent power laws are similar to those 
found for compact patterns. The density of cancer cells 
and their average division rates through the growth pat- 
terns are not significatively altered when the tumors are 
submitted to mild treatments (long time intervals be- 
tween consecutive drug doses). Moreover, their growth 
patterns are very similar to those found in the untreated 
counterparts. Finally, tumors exhibiting necrotic cores 
are more easily eliminated when shorter periods of the 
drug administration are considered (Fig. 0). 




x=l 



T-2 



T=4 



1=5 



t=10 



t=20 



FIG. 5: Compact patterns under treatment with antimitotic 
drugs. The dose period r was varied and 9^ =0.3 was 
fixed. The other parameters are indicated in table |. Tumors 
become more and more fractal as the dose period decreases. 
At a critical period the tumors reach a frozen state in which 
their size remain constant during all the treatment. 



Several drugs used in cancer chemotherapies do not kill 
cancer cells. Instead, they aim to stop the cell cycle in 
specific checkpoints. As a consequence of that, the tu- 
mors cease or slow down their growth rate. Examples of 
such drugs include the antimitotic agent Curacin A that 
blocks the cell cycle in the M-phase, mitomycin C, dox- 
orubicin, and aclarubicin, among others [|29|. In order 
to analyze the effects of drugs that inhibit cell division 
on tumor patterns, we introduce a very simple change in 
the cell dynamics of the model described in the previous 
section. Instead of an additional death probability given 
by Eq. (^) , the division probability is modified by includ- 
ing an additional term in Eq. (^). So, the new division 
probability becomes 



Pdiv = 1 - exp 



N 



Q 



9 (Q) 

div , 



(10) 



in which the first and second terms in the exponential 
argument compete between themselves stimulating and 
inhibiting the cell mitosis, respectively. The parameter 
^ctiv con trols the drug influence on the cell division. Ob- 
viously, if the exponential argument is larger than the 
unity then Pdiv = 0. 

The biological interpretation of the present model is 
very distinct from that considered in the last section. 
Now, the drugs modify the intrinsic characteristics of the 
cancer cells and, consequently, the tumors must behave 
differently from those treated with cytotoxic agents. In- 
deed, as one can see in Fig. ^ a morphological transition 
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4 8 12 16 20 

X 



FIG. 6: Critical values 9 as a function of the dose period r. 
The symbols represent simulational data and the straight lines 
the respective exponential fittings. The open circles refer to 
the simulations using the parameters for the compact patterns 
listed in table Q. The squares (filled circles) represent the 
same simulations except that T x — 6 (&div = 0.1). Also, the 
slopes 1/4 (circle) and 1/6 (square) are indicated. In order 
to estimate we realized 20 independent simulations for 

each value of 0^9' and considered that the growth failure for 
all these tentatives means tumor latency. 



for the tumor patterns when the dose periods decrease is 
observed . If the period is sufficiently short the tumor size 
remains constant during the therapy. Obviously, the crit- 
ical period depend on other model parameters, specially 
on 6* £ . In order to see this dependence, was evaluated 

the critical value of 6^* for which the tumors cease their 
growth with probability 1 as a function of t. A mean 
field analysis of Eq. (|l(J) provides 6^ ~ exp(— r/T x ), 
in agreement with the simulations as indicated in Fig. 0. 
Thus, it was found an exponential decay with an univer- 
sal characteristic length T x that is independent of the 
other model parameters. This law, relating two impor- 
tant clinical parameters, is valid for an wide set of param- 
eter values. In addition to the compact tumors shown in 
figure ||, it was also simulated papillary tumors under an- 
timitotic treatment, and a similar morphology transition 
was observed. 

Fig. [t] shows that these morphological transitions oc- 
curred at well defined offl values. Near the critical value, 
instabilities in the reaction-diffusion equations lead to 
a branching growth Below this threshold the cell 

death rate equals cell division, and the tumor growth ef- 
fectively stops. This qualitative behavior near the tran- 
sition will probably not be observed in experimental tu- 
mors, since a very reliable determination of 9\V is very 
difficult. On the other hand, experiments exhibiting both 
regimes, above and below the critical value, are realiz- 
able today, due to the great advances in drug delivery 
control. Such assays might corroborate the complex be- 
haviors predicted by our model. 

The antimitotic treatment affects both the growth and 



41 * 

0.25 0.26 0.28 0.29 

FIG. 7: Morphological transitions at . The patterns 

were originally compact and the model parameters are listed 
in table The dose period was fixed at r = 2 and 6^ varied 
around its critical value. 



scaling laws. The compact patterns in the original model, 
for which the gyration radius and number of occupied 
sites on the tumor border scale as the square root of the 
total number of occupied sites (i. e. v = a = 1/2), be- 
come progressively more fractal as the dose intervals are 
decreased. Consequently, it was obtained v > 1/2 and 
a > 1/2 indicating fractal tumors with a fractal dimen- 
sion given by df — 1 jv. In general, all the tumors become 
more fractal when submitted to antimitotic treatments. 
So, there is a neat contrast with the invariance of the 
fractal dimension of the growing tumors under the cy- 
totoxic treatment described in the previous section. In 
Fig. H the increase in the number of cancer cells relative 
to the treatment beginning (Nc — No) is plotted, for dis- 
tinct dose intervals, as a function of the time afterwards 
the treatment beginning (T — To). These curves suggest 
power laws for the time evolution of the number of cancer 
cells with a weak dependence on the r parameter. Indeed, 
we have (N c - N ) ~ (T — T o ) whit (3 G [1.2, 1.4]. As 
r decreases, the number of tumor cells initially decays 
for a while, but subsequently recovers its growth (inset 
of Fig. |[). Below a given dose interval, Nc decreases 
monotonically up to a constant value. 

In order to analyze the cell division rates through 
the tumor, we computed the average cancer cell density 
((<7 C )) and cell division rate (wdiv) along a longitudinal 
cut across the growth pattern. These average values were 
plotted as a function of the distance from the capillary 
vessel. In Fig. || these plots are shown for compact pat- 
terns treated with distinct r values. For larger r values 
(Fig. §(a)) both, division rate and cancer cell density, 
have sharp maxima on the tumor border in front and 
opposite to the capillary vessel. Notice that the peaks 
for division rate are significatively narrower than those 
for cancer cell density, demonstrating that the prolifer- 
ative fraction of the tumor comprises just a small rim 
located at the tumor border. As the treatment is inten- 
sified (shorter r values are considered) the cell densities 
through the tumors become more uniform, whereas the 
division rates continue exhibiting sharp maxima at the 
tumor borders (Figs, ^(b) and (c)). If r is sufficiently 
short in order to halt the tumor growth, the division 
rates and cancer densities become uniform along the tu- 
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FIG. 8: Growth curves (dashed lines) for the papillary pat- 
terns under antimitotic therapy. The dose intervals used were 
r = 20, 10, 7, 6, 5, 4, 3, 2 and the upper curves represent the 
larger intervals. For the larger dose intervals (r > 5) the tu- 
mor growth follows power laws with exponents in the range 
[1.2,1.4]. For shorter intervals (r < 4), the curves are shown 
in a semi-log plot (inset). In these simulations 6^ = 0.3 was 
used in addition to the parameters of table |. 



mor (Fig. |](d)). So, in these cases the division rates are 
counterbalanced by the death rates, leading to a vanish- 
ing net rate of cancer growth. All these results show that 
the tumor patterns are altered when the mitotic proper- 
ties of the cancer cells are modified by external agents. 



V. CONCLUSIONS 

In the present work, a reaction-diffusion model to sim- 
ulate the effects of chemotherapies on the growth of carci- 
nomas "in situ" have been studied. The model includes 
cell death and division, competition among cancer and 
normal cells by nutrients and periodical drug adminis- 
tration. Two kinds of chemotherapies, using cytotoxic 
and antimitotic drugs, were modelled. 

In the cytotoxic model the tumors can be completely 
eradicated, cease their growth or grow continuously. The 
last case occurs when the treatment is inefficient or the 
intervals between consecutive doses of drugs are large. 
Moreover, the morphologies and scaling laws of the grow- 
ing tumor patterns are preserved. In contrast, for ther- 
apeutic approaches using antimitotic agents a morpho- 
logical transition in the tumor patterns was observed. 
The growth patterns become progressively more fractal 



as more effective treatments (shorter intervals between 
consecutive doses and/or more efficient drugs) are con- 
sidered. 

Such morphological transitions are similar to re- 
cent studies claiming that bacterial colonies exposed to 
non-lethal concentrations of antibiotics exhibit drastic 
changes in their growth patterns |3^| . For bacteria, these 
changes were imputed to variations in bacterial proper- 
ties such as metabolic load and chemotaxis. In turn, 
normal and cancer cells cultured in monolayer and colla- 
gen gel exhibit a dynamical transition in their aggrega- 
tion regimes as an adaptive response to the growth con- 
straints imposed by high cell population density or long 
permanence in culture |33 . Again, the results obtained 



in the present work are in agreement with the point of 
view that cancer cells can develop an integrated defense 
program against stress situations similar to the response 
of bacterial colonies facing severe and sustained threats 
|nj. However, as far as we are concerned, there are no 
reports on cancer literature concerning transitions mor- 
phologies in the histological patterns of tumors submitted 
to chemotherapy. Thus, formal models like the one pro- 
posed in this paper, generally not familiar for the most 
of the biomedical researchers |34| , can guide and refine 
new experiments intended to analyze such morphological 
transitions. Currently, experiments "in vitro" with can- 
cer cells are being performed in our laboratories in order 
to investigated those morphological transitions. 

It is important to mention that, usually, real effective 
treatments employ several therapeutic methods simul- 
taneously [^Tj. So, a more realistic chemotherapeutic 
model should consider combined cytotoxic and antimi- 
totic treatments. Nevertheless, in these cases relevant 
features such as the robustness of the tumors fractal scal- 
ing and the morphology transitions can be masked. At 
last, we are modelling others cancer therapeutic strate- 
gies, such as virus and immunotoxins therapies, as well as 
the combination of distinct treatments by using reaction- 
diffusion models like that considered in this paper. 
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